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We estimate the production of two meson pairs in high energy photon photon collisions produced 
in ultraperipheral collisions at LHC. We show that the study of charge asymmetries may reveal the 
existence of the perturbative Odderon and discuss the concrete event rates expected at the LHC. 
Sizable rates and asymmetries are expected in the case of proton-proton collisions and medium 
values of 77 energies ^/s 77 w 20 GeV. Proton-proton collisions will benefit from a high rate due to a 
large effective 77 luminosity and ion-ion collisions with a somewhat lower rate from the possibility 
to trigger on ultraperipheral collisions and a reduced background from strong interactions. 
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I. INTRODUCTION 

Hadronic reactions at low momentum transfer and high energies are described in the framework of QCD in terms of 
the dominance of color singlet exchanges corresponding to a few reggeized gluons. The charge conjugation even sector 
of the t— channel exchanges is understood as the QCD-Pomeron. The charge-odd exchange is less well understood, 
although the need for the Odderon contribution [H, @j j m particular to understand the different behaviors of pp and pp 
elastic cross sections Q, is quite generally accepted. Studies of specific channels where the Odderon contribution is 
expected to be singled out have turned out to be very disappointing 4] but new channels have recently been proposed 
[ETS 0, Hi ■ Another strategy to reveal the Odderon has been first initiated in Ref. @ , where it was stressed that the 
study of observables where Odderon effects are present at the amplitude level - and not at the squared amplitude level 
- is mandatory to get a convenient sensitivity to a rather small normalization of this contribution. This approach has 
been extended to the production of pion pairs [TH EH EH EH EH EH since the ir + n~ -state does not have any definite 
charge parity and therefore both Pomeron and Odderon exchanges may contribute to the production amplitude. 

In this paper, we study within perturbative QCD (pQCD) the charge asymmetries in the production of two pion 
pairs in photon-photon collisions 

7 (<z,e) 1 (q',e')^ir+(p + ) 7r+(p' + ) (1) 

where e and e' are the initial photon polarization vectors, see Fig. [TJ We have in mind the ultraperipheral collisions 
of protons or nuclei of high energies, which constitute a promising new way to study QCD processes initiated by 
two quasi real photons [la, EH EH- At high energy the application of pQCD for the calculation of a part of this 
process is justified by the presence of a hard scale: the momentum transfer t = (q — p + — p-) 2 from an initial photon 
to a final pion pair, —t being of the order of a few GeV 2 . The amplitude of this process may be calculated within 
fey-factorization, as the convolution of two impact factors representing the photon to pion pair transitions and a two 
or three gluon exchange. The impact factors themselves include the convolution of a perturbatively calculable hard 
part with a non-perturbative input, the two pion generalized distribution amplitudes (GDA) which parametrize the 
quark- antiquark to hadron transition. Since the 7r + 7r~ system is not a charge parity eigenstate, the GDA includes 
two charge parity components and allows for a study of the corresponding interference term. The relevant GDA is 
here just given by the light cone wave function of the two pion system fl9j . 

II. KINEMATICS 

Let us first specify the kinematics of the process under study, namely the photon-photon scattering as stated in 
Eq. |T]). We introduce a Sudakov representation of all particle momenta using the Sudakov light-like momenta pi,p2- 
The initial photon momenta are written as 
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Figure 1: Kinematics of the reaction 77 — > tv + tv n + n in a sample Feynman diagram of the two gluon exchange process. 



where s = 2p\ ■ p 2 . The momenta of the two pion systems are given by 

u _ /, Pl-K v ji , m L + & a , u 2 _ _ -2 

P2?r — U ~/>Pl H I ^2 +^27r±' ^27r± — ^27r ; 



(3) 



The (massless) quark momentum ?! and antiquark momentum l 2 inside the upper loop before the formation of the 
two pion system are parametrized as 



iu u 1 (J 1 z P2tt) u n \u 
I'l = Zpl H p% + {Ix + ZfrsnXr, 



(4) 



1% = zpl 



(-l + zp 2v ) 2 M 



Pi + ( _ ^-L + ZP2nxT, 



(5) 



where 21 is the relative transverse momentum of the quarks forming the two pion system and z — 1 — z, up to small 
corrections of the order p 2n /s. Following the collinear approximation of the factorization procedure in the description 
of the two pion formation through the generalized distribution amplitude, we put / = in the hard amplitude. 
In a similar way as in ^ , ^ we parametrize the momenta of the produced pions as 



/( , ,x . ml + {p + CP2^) 2 p, / , t ,„ 

p+ = Cpi + - r P2 + (p± + (P2ir±r, 



(6) 



P- = (Pi 



ml + {-p + (p2ir)' 



P 2 l + (-p± + (P2*±T, 



(7) 



where 2p is now their relative transverse momentum, £ 



is the fraction of the longitudinal momentum p 27r 



carried by the produced 7r + , and £ = 1 — £. The variable £ is related to the polar angle 6 which is defined in the rest 
frame of the pion pair by 



/3cos0 = 2C-l, 



4 ml 



(8) 



Since the "longitudinal part" of the two pion wave function depends only on the angle and does not depend on 
the azimuthal decay angle <f> (in the same rest frame of the pair), we focus on the calculation of charge asymmetries 
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expressed in terms of 9 (see below). Similar expressions as Eqs. (|4][8]) are used for the lower quark loop. Since we 
are interested in photon interactions in ultrapcriphcral collisions of hadrons, the photons are quasi real and hence 
predominantly transversely polarized. The polarization vectors of the photons are written as : 

e(+) = -^=(l, l ), %-) = ±(l,-i). (9) 

We will consider spin averaged cross sections since hadron colliders do not produce polarized photon beams. 

III. SCATTERING AMPLITUDES 

It is well known that at high energies (s ^> \t\) the amplitude factorizes into impact factors convoluted over the 
two-dimensional transverse momenta of the i-channel gluons. 

For the Pomeron exchange, which corresponds in the Born approximation of QCD to the exchange of two gluons 
in a color singlet state, see Fig. [TJ the impact representation has the form 

M r = -is / \2Z2Z2 J p {ki,k 2 )- j; r (fci,fc 2 ), (10) 

where Jp T (fci, k 2 ) is the impact factor for the transition 77- — * 7r + 7r~ via Pomeron exchange. 

The corresponding representation for the Odderon exchange, i.e. the exchange of three gluons in a color singlet 
state, is given by the formula 

8tt 2 s f d% d 2 k 2 d 2 k 3 S^jh + k 2 + k 3 - p 27r ) r -. -. -> -. 

Mo = -— / ^t2tit2 ^> (fci , k 2 , k 3 ) ■ (ki , k 2 , k 3 ), (11) 

J ^Z7TJ rZ<2 ft 3 

where Jq t {k\, k 2l k 3 ) is the impact factor for the transition 7^ — > tt + tt~ via Odderon exchange. 

The impact factors are calculated by the use of standard methods, see e.g. Ref. [2(| and references therein. 

A. Impact factors for 7 — > -k + -k~ 

The leading order calculation in pQCD of the impact factors gives 

2 b 1 

J^ T (h,k 2 ) = - ie A 9 J a f dziz-^^-Qrikuk^^izX,™ 2 ^), (12) 



4N C 







where the vector Qp(ki, k 2 ) is defined by 



- ■ ,- zp 27r zp 2lT kx~zp 2v h-zp 2v 

Mki, k 2 ) = ^ — + — — — — — -. (1A) 

z 2 pl + ^ 2 z 2 P 2 \ + V 2 {h- zp^y + fi 2 {ki-zp 2 ^) 2 +n 2 



The calculation of the Odderon exchange contribution gives 



1 

i p n 3 r\ abc f -> -. _ - 1 



JS r (k 1 ,k 2l k 3 ) = / dz(z-z) e{T) -QoikxMM -$ i=0 (z,C,m^), (14) 





where we have used the definition 



t(ki,K 2 ,k 3 ) = -5-^5 ^ + 37^ o + > — h — ■ (15) 

z 2 pl + v 2 z 2 pI + ^ 2 ^\( ki -zp 27r ) 2 + fi 2 {ki-zptoY + f) 
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B. Generalized two pion distribution amplitudes 



A crucial point of the present study is the choice of an appropriate two pion distribution amplitude (GDA) [12], Il3l . 
ITU . liH [H, H3, H3, Hl| which includes the full strong interaction related to the production of the two pion system. 

Based on an expansion of the GDA in Gegenbauer polynomials C™(2z — 1) and in Legendre polynomials P;(2£ — 1) 
[24l . l2l| it is believed that only the first terms give a significant contribution: 



$ 7 1 {z,(,m 2n ) = 6zz/3/i(m 27r ) cos6», (16) 



$ /_0 (z, C, m 2 Tr) = 5zz(z-z) 



3 - 8 2 

^/o(m 27r ) + /3 2 /2(m 2T )P 2 (cos0) 



(17) 



where /i(m 27r ) can be identified with the electromagnetic pion form factor i r ' 7r (m 2w ). For our calculation we use the 
following i^-parametrization inspired by Ref. [26] 

AK) = fxK) = ^^(O ^yi (is) 



with 



to* 



fl^K) = -^—2 J~~ F 7-2 ; (19) 



2/2 x 9 \ 3/2 

o ; _ p '"„ f - Ami X 
TO? — 4m 2 , 



2 

As masses and widths we use m p = 775.49 MeV, T p = 146.2 MeV, m w = 782.65 MeV, I\, = 8.49 MeV [27|. We fit the 
remaining free parameters to the data compiled in Ref. [28| obtaining a = 1.78 • 10~ 3 and b = —0.154. Including a 
hypothetical p' resonance as originally used in Ref. [2(| , gives a significant better fit to the data at large m 2Tr but has 
only small effect on the asymmetry which is the main object of our studies. 

For the 1 = component we use different models. The first model follows Ref. [l2j and reads 

/o/ 2 (m 2w ) = e ^(™-) \BW h/ Ml«)\ ■ (22) 

The phase shifts 5 / 2 are those from the elastic tt + tt~ scattering, for which we use the parametrization of Ref. [5^] 
below and that of Ref. jlc| above the KK threshold. 

\BWf 0/2 (m 27V )\ is the modulus of the Breit-Wigner amplitudes 

2 

BWf 0/2 (m 2 2w ) = — 2 , (23) 



with m/o = 980 MeV, r /o = 40- 100 MeV, to^ = 1275.1 MeV, T h = 185 MeV [27 1. 

In the second model - elaborated in Ref. 23j - the functions / / 2 are the corresponding Omnes functions for S— 
and D— waves constructed by dispersion relations from the phase shifts of the elastic pion scattering: 

Mm*,) = exp ( ttJ, + ^ r ds 2( Sl{s 2 ] , ) , with k = i ^ ds^. (24) 

The assumption that the phases of the GDA equal those of the elastic scattering looses its solid base beyond the 
KK threshold. As discussed in Refs. [H|, [3l| it might well be that the actual phases of the GDA are closer to the 

2i5 l 

phases of the corresponding T matrix elements m e 2i . The third model for the 1 = component of the GDA takes 

this into account by using the technique of model 2 with these phases 5t,i of the T matrix elements. 

While the first and the second model give quite compatible results, model three differs from them significantly. 
The most striking difference is the absence of pronounced fo resonance effects in model 3. In fact, measurements 
at HERMES [12] do not observe a resonance effect at the /o-mass even though a confirmation by an independent 
experiment would be desirable. From the same measurements at HERMES one can draw the conclusion that using 
St,2 and 5 2 for the f 2 region are both compatible with data 23]. Having this in mind, we consider also a fourth model 
- a mixed description with the fo contribution from model 3 and the f 2 contribution from model 2. 
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C. Photon exchange amplitude 

The photon has the same C-parity as the Odderon. Therefore, its exchange between the two quark-antiquark 
systems can mimic an Odderon exchange. The according amplitude is straightforward to calculate and reads 

M^ = ±J?'J?, (25) 

with 

J ? = i j ^(*- 2 )^-^(^+^ + j?+W)* z ^ ,c ' <) - (26) 



In our concrete process the photon exchange amplitude amounts to the order of 10% of the Odderon exchange 
amplitude. Although we do not neglect this contribution, it is clear that the asymmetry described in the following 
section is driven by the Odderon/ Pomcron-interference. 



IV. CHARGE ASYMMETRIES AND RATES 



In the following we make use of the fact that the GDAs for C-even and for C-odd pion pairs are orthogonal to each 
other in the space of Legendre polynomials in cos 9. As a consequence, in the total cross section the interference term 
completely vanishes. In contrast only the interference term survives, when the amplitude squared is multiplied by cos 9 
before the angular integration which corresponds to selecting the charge asymmetric contribution. The asymmetry 
we are interested in is defined as 



A(t,ml„,m^) 



j cos 9 cos 9' da(t, m^., m^, 6>, 6*') 

(27) 



/ da(t,m 2 2n ,m' 2 \,9,9') 
/_ x dcos6> J_ x dcos9' 2cos# cos9' Re [Mp{M® + M-,)* 



d cos 9f\d cos 6" \M P \ 2 + \M® + M 



■ I 2 



The obtained landscape as a function of the two invariant masses is not particularly illuminating and would be 
difficult to measure. To reduce the complexity, we integrate over the invariant mass of one of the two pion systems 
to obtain 



A{t, m^Tr'i TO min' TO max) 



dm% J cos 9 cos 9' da(t, m%„, m'^,9, 9') 
Chdrngf da(t, m t„, m g,9,9>) 

dmgj\dcos9 J\d cos 9' 2cos<9 cos 9' Re [M V {M® + M y )*\ 



f™ 2 — dm'l dcos9f* 1 dcos9' \M P \ 2 + \M® + M 



. I 2 

'7 



(28) 



Let us note that since two pion pairs are always in the same C-parity state, because of C = + parity of the 
initial 77 state, it is necessary to keep in Eq. (f55|) the integration weight cos6>'. Without this weight the single 
charge asymmetry would vanish. The deviations from this vanishing of the asymmetry may serve as a measure of 
experimental uncertainties. 

Considering an analytic calculation of the matrix element in Eq. (|lip . it turns out that it would require the 
calculation of two dimensional two-loop box diagrams, whose off-shellness for all external legs is different. The 
techniques developed in Refs. [HI, [34], [3]| can not be applied here due to the very elaborated topology of the most 
complicated diagrams involved (square box with one additional diagonal line). Instead we rely on a numerical 
evaluation by Monte Carlo methods. In particular we make use of a modified version of Vegas as it is provided by 
the Cuba library (36|. The result for the asymmetry A at t = — 1 GeV 2 (resp. t = —2 GeV 2 ) is shown in Fig. [2] (resp. 
Fig. [3]). Since our framework is only justified for m| ff < —t, (in fact strictly speaking, one even needs m| ff <§C —t ), we 
keep 7712,7 below 1 GeV (resp 1.4 GeV). In each case, we present the expected asymmetry with the GDAs parametrized 
as discussed above. 

In order to evaluate the feasibility of the Odderon search, we need to supplement the calculation of the asymmetry 
with rate estimates in ultraperipheral collisions at hadron colliders. This rate depends on the Pomeron dominated 
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Figure 2: Asymmetry A at t = — 1 GeV 2 for model 1 (solid), 2 (dashed), 3 (dotted), and 4 (dash-dotted) - model 3 and 4 
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Figure 3: Asymmetry A at t = —2 GeV 2 for model 1 (solid), 2 (dashed), 3 (dotted), and 4 (dash-dotted). Left column has 
m m i n = .3 GeV and m max = m p , while right column has m m i n = m/ and m ma x = 1.4 GeV. 



photon-photon cross section and on the equivalent photon flux. The total photon-photon cross section falls off rapidly 
with increasing \t\ (see Fig. 2]). Therefore, the integration mainly depends on the lower limit of |i|-integration (i m ; n ). 
Already for t min = — 1 GeV 2 we find er 77 = 1.1 pb. 

Although in Ref. |16| it is claimed that the photon flux is best for medium-weight ions, and especially superior to 
that of protons, this is in fact not true. To obtain Fig. 6 of Ref. [l6| the effective 77 luminosity has been calculated 
for protons and ions by the Monte Carlo program Tphic [37J| which is based on the Weizsacker- Williams method 
[38l l39j] with the additional condition of non-overlapping ions [4Cl l4ll] . but the authors used a quite small luminosity 
for proton-proton collisions at the LHC (14 000 mb _1 s _1 instead of the official design luminosity 10 7 mb _1 s^ 1 [42l. |43| ) . 
Very unfortunately, the identical figure is reprinted in Ref. (l7j while a non consistent p-p luminosity of 10 7 mb _1 s _1 
is cited. 

As already was shown by Cahn and Jackson [4lj . the 77 luminosity in case of ions can be expressed in terms of 
a universal function £(z), where z = MR « 5.665A 4 / 3 with M being the mass and R the radius of the ion, and a 
prefactor proportional to Z 4 . Since the luminosity for ions at the LHC decreases roughly as Z~ 4 , the prefactor's Z 
dependence is more or less compensated and only the universal function £ remains which is exponentially decreasing 
with z. Hence, lighter projectiles provide a larger effective 77 luminosity, with the protons offering the highest 
luminosity. 

Therefore, we provide a corrected overview over the various effective 77 luminosities in Fig. [5J For the different ion 
scenarios that are discussed in Ref. [44[ we use the parametrization of Ref. l4ll which relies on the Weizsacker- Williams 
method [H,[3^] with the additional condition of non-overlapping ions [40ll4lj. For protons usually a calculation based 

on the proton electric dipole form factor Fe(Q 2 ) = 1/(1 + Q 7 ® Gc y-i ) m combination with the Weizsacker- Williams 
method is used, as it can be found in Ref. [45j . In Ref. [46| a slightly imp roved version is given, which lowers the 
photon flux. An inclusion of the corresponding magnetic dipole moment [47] would lead to a flux between those both. 
For this reason, we use the formulas given in Ref. [45| and Ref. [!(| for the case of proton-proton collisions. We 
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Figure 4: t-dependence of photon-photon cross section. 
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Figure 5: Effective 77 luminosities for the collision of p-p based on Ref. [45| (dash-dotted) and Ref. [4(| (dashed). The results 
using the parametrization of Ref. [4l| for ions are given by solid lines for p-p, Of 
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Pblos f rom t°P to bottom. For ions we used the average luminosities as given in Ref. [44] , for proton we used L pp — 10 34 cm z s 



For comparison also effective 77 luminosities at the ILC are given for ^/s e + e - = 250 GeV and y/s e + e 
dotted lines). 



500 GeV (both as 



also provide the results for the proton treated as a heavy ion because it might be that the non-overlap condition - 
reducing the flux - is of importance [4q , even if such a procedure does not include the proton form factor. We consider 
such a lower result as a conservative lower estimate. Since the luminosity factors will cancel in the asymmetry, this 
uncertainty does not affect our conclusions on the Odderon effects. For comparison, we provide also the effective 77 
luminosities at the intended ILC where the design luminosity for e + e~ collisions is 2 • 10 34 cm _2 s _1 [48| 

As shown in Fig. the effective 77 luminosity decreases rapidly with increasing energy. Since our intermediate 
hard scale t is quite small, we are not forced to consider extremely large photon-photon energies. The only condition 
on the minimal photon-photon energy (s m ; n ) is s mm \t] to ensure the validity of high energy factorization into two 
separated pion systems. Hence, a s mnl of 400 GeV 2 could be already enough. The effect of varying s m ; n is displayed 
in Fig.© 

For p-p collisions, the rates are high and even for < m ; n = —2 GeV 2 sizable. Although heavy ions would offer the 
possibility to trigger on ultraperipheral collisions by detecting neutrons from giant dipole resonances (GDR) in the 
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100 GeV. 
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Figure 6: Rate of production of two pion pairs in ultraperipheral collisions in dependence on the lower cut s m i n given in 'events 
per month' in case of ions, and 'events per six months' in case of protons which in both cases correspond to one year of running 
of LHC. The l.h.s. plot shows the rate for t m i n = — lGeV 2 , the r.h.s. that for for t m i n = — 2GeV 2 . The solid line displays 
the result for p-p collision using Ref. [HI, the dashed-dotted that for protons treated as heavy ions, the dashed one that for 
Ar-Ar collisions, and the dotted line that for Pb-Pb collisions. On the left figure, also the much smaller rates coming from the 
Odderon exchange are shown (with the same dashing). 

Zero Degree Calorimeters, the rates that can be read off from Fig. [6] are rather low. Only for medium-weight ions 
there might be the possibility to measure the process. 

In hadron-hadron collisions the process of interest could as well be connected by Pomerons to the colliding hadrons. 
In that case one would have to deal with the unknown hadron-Pomeron and Pomeron-Pomeron-two-pion couplings. 
Therefore, we consider that process as background. This circumstance would make heavy ions preferable because 
of the different scaling of Pomeron and photon coupling when changing from proton to nuclei scattering. However, 
experimentally such a background can be suppressed by refusing events with a total px larger than some small cut-off. 
Indeed, the 77 events dominates, due to the photon propagator singularities, when each of the transverse momentum 
is small, which on the average is satisfied when the total transverse momentum is imposed to be small. 

Inclusion of high energy evolution of the two gluon exchange a la BFKL [5(1 HH [H, HE HH [Efl would increase the 
total cross section but at the same time diminish the asymmetry since the Odderon including BKP evolution [56|, [57| 
has a smaller intercept. The study of these effects goes beyond the goal of this paper. Since we demonstrated that 
the rate and asymmetry were sizable for moderate values of the 77 energy, we do not expect evolution effects to play 
a major role in the proposed strategy to uncover Odderon effects. 

V. CONCLUSION 

We have investigated in real photon-photon collision the production of two 7r + 7r~-pairs well separated in rapidity. 
Due to the non fixed C-parity of these pairs, beside a Pomeron exchange an Odderon exchange is possible as well. We 
have calculated both contributions in a perturbative approach - justified by t providing the hard scale - where the 
only soft building blocks needed are the GDAs of the pion pairs. We have shown that a charge asymmetry in the polar 
angle 9 (defined in the rest frame of the pion pairs) is linearly dependent on the Odderon amplitude and moreover 
is sizable but GDA-model dependent. In fact, the predicted asymmetry depends much of the two pion GDAs, which 
are still not really known although HERMES measurements of two pion electroproduction [32j disfavor models with 
a strong /o coupling to the 7r + 7r~ state. We however think that higher statistics data, which may come from a JLab 
experiment at 6 or 12 GeV, are definitely needed before one can trust a definite model of the GDAs. Because two 
pion deep electroproduction in the low energy domain is dominated by quark exchanges, this test of the GDA models 
is independent of any Odderon search. This looks like a prerequisite to a trustable extraction of the Odderon signal 
- or of an upper bound on Odderon exchange amplitudes - from ultraperipheral collisions. 

We have discussed the possibility to measure this process in ultraperipheral collisions at LHC. Although in ion- 
ion collision it is easier to trigger on ultraperipheral collision with a lower background from diffractive events from 
strong interactions, the event rates are too low. In proton-proton collision the rates are high enough to measure 
the asymmetry even though isolating ultraperipheral collisions in proton-proton collisions may be a challenge for 
experimentalists. 

Let us finally note that the background from strong interactions would be completely absent in an e + e~ collider 
which via Compton-back-scattering could work as a very effective 77 collider. At the ILC [48| for a nominal electron 
beam energy of 250 GeV the luminosity for photon-photon collisions would be 3.5 • 10 33 cm -2 s -1 [58[ for the high 
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energy photons (energy fraction at least 80% of the maximal possible energy fraction) or even higher, if optimized for 
photon-photon collisions [591 ] - As Fig. O reveals, even running in the electron-positron mode the effective 77 luminosity 
is slightly larger than at the LHC. For these reasons, the ILC would provide an ideal environment to study the process 
of interest. 
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